home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 005 / miscstat.arc / FDIST.BAS < prev    next >
Encoding:
BASIC Source File  |  1986-04-25  |  1.1 KB  |  43 lines
  1. 1 'F-DISTribution
  2. 2 KEY OFF:SCREEN 0,1,0:WIDTH 40:COLOR 15,1,9:CLS
  3. 10  PRINT "F-DISTRIBUTION"
  4. 20  PRINT 
  5. 30  PRINT "(TO END PROGRAM ENTER AN F-VALUE OF 0)"
  6. 40  PRINT "F-VALUE";
  7. 50  INPUT F:IF F<0 THEN 50
  8. 60  IF F = 0 THEN 340
  9. 70  PRINT "DEGREES OF FREEDOM IN NUMERATOR";
  10. 80  INPUT D1
  11. 90  PRINT "DEGREES OF FREEDOM IN DENOMINATOR";
  12. 100  INPUT D2
  13. 110 X = 1
  14. 120  IF F < 1 THEN 170
  15. 130 S = D1
  16. 140 T = D2
  17. 150 Z = F
  18. 160  GOTO 200
  19. 170 S = D2
  20. 180 T = D1
  21. 190 Z = 1 / F
  22. 200 J = 2 / 9 / S
  23. 210 K = 2 / 9 / T
  24. 220 Y =  ABS ((1 - K) * Z ^ (1 / 3) - 1 + J) /  SQR (K * Z ^ (2 / 3) + J)
  25. 230  IF T < 4 THEN 270
  26. 240 X = .5 / (1 + Y * (.196854 + Y * (.115194 + Y * (.000344 + Y * .019527)))) ^ 4
  27. 250 X =  INT (X * 10000 + .5) / 10000
  28. 260  GOTO 290
  29. 270 Y = Y * (1 + .08 * Y ^ 4 / T ^ 3)
  30. 280  GOTO 240
  31. 290  IF F >  = 1 THEN 310
  32. 300 X = 1 - X
  33. 310  PRINT "PERCENTILE: ";1 - X
  34. 320  PRINT "TAIL-END VALUE: "X
  35. 325  PRINT 
  36. 330  GOTO 40
  37. 340 PRINT:LOCATE 24,1:PRINT"PRESS <SPACE> TO RUN AGAIN";
  38. 350 LOCATE 25,1:PRINT"PRESS <ESC> TO EXIT";
  39. 360 K$=INKEY$:IF K$="" THEN 360
  40. 370 IF K$=" " THEN RUN
  41. 380 IF K$<>CHR$(27) THEN 360
  42. 10000 RUN"MENU
  43.